ASPECTS OF STABLE POLYNOMIALS 



STEVE FISK 



This article is an introduction to the properties of stable polynomials in 
several variables with real or complex coefficients. These polynomials are 
defined in terms of where the polynomial is non-vanishing. 

'All polynomials /(xi, . . . , x^) with complex coeffi- 
cients such that /((Ti, . . . , ad) 7^ for all ai, . . . , ad 
in the right half plane. If we don't need to specify d 
,we simply write 7^(C). 
We call such polynomials stable polynomials. In one variable they are often 
called Hurwitz stable. 



Definition 1. HdiC) 



This article is a companion to [7] where we studied polynomials non- 
vanishing in the upper half plane. 

'All polynomials f{xi,...,Xd) with complex coeffi- 
cients such that /(cJi, . . . , ad) / for all cji, . . . , aa 
in the upper half plane. If we don't need to specify d 
,we simply write U (C). 
Such polynomials are called upper polynomials. 



Definition 2. (C) 



Since multiplication by i takes the right half plane to the upper half plane, 
we have the important fact that 

/(xi, ...,Xd)e\]d (C) <^ f{txi, txd) G HdiC) 

Consequently, many of the properties of 7i{C) follow immediately from the 
properties of U (C), and their proofs will be omitted. See [THTlfTS] for more 
information about U (C) and 7i{C). We use the notation in [3, and let RHP 
denote the right half plane. 

We do not cover well-known topics in one variable such as Routh-Hurwitz, 
the Edge theorem, and Kharitonov theory. See |lHll2j. 



1. Complex coefficients 

These results are proved following the corresponding arguments for upper 
polynomials. The main difference is the disappearance of minus signs for 

n{c). 

Fact 1. 

(1) Suppose /(x) G HdiC). 

(a) Ifa^O then af{x) G Hd(C). 
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(b) If ai > 0,. . . ,ad > then /(aixi, . . . , a^Xd) € HaiC). 

(c) // 3?(c7i) > 0, . . . , 3f?(crd) > then f{xi + ai, . . . ,Xd + aa) G 

(d) //3?(ct) >0 then f{a,X2,...,Xd) G W(C). 

(e) /(xi + y, X2, . . . , Xrf) G W(C). 

(f) /(x,x,X3, . . . ,xd) G 7Y(C). 

(2) // /(x) G n{C) then f{ta, X2, . . . , x^) G H(C) U {0} i/ a G M. 

(3) H (C) is closed under multiplication and extracting factors. 

(4) IfEoM^)y' e H(C) i/ien Eo /^Wy""' G W(C). 

(5) ///(x) G n{C) then a,J(x) G H(C) U {0}. 

(6) IfT.fi{^)y' e n{C) then all coefficients fi{x) are in 7^(C) U {0}. 

We have interlacing for both 7i(C) and U (C); later we will see there is 
another kind of interlacing when all the coefficients are real. Note that 
interlacing is symmetric for stable polynomials. 

Definition 3. 

/(x) ^ <7(x) if and only if /(x)+y<7(x)G n{C) 
/(x) ^ 5(x) if and only if f{^)+yg{^)e U (C) 

Fact 2. 

(1) If Eo n(.^)y' e H{C) then /,(x) ^ /,+i(x) /or i = 0, . . . , n - 1, 

provided fi and /j+i are not both zero. 

(2) ///GH(C) thenfJ^d^Ji^). 

(3) Suppose f,gen{C) 

(a) fgJLghiffhe n{C) and f g. 

(b) IffJLgtheng^f. 

(c) 7/ / < — 5 anc? / < — /i i/ien / < — g + h. 

(d) // / and /i g then f + h 5. 

(e) If f ^ g ^ h then f + h^ g. 

(4) (Hermite-Biehler) Suppose that /(x) is a polynomial, and write /(x) = 
/e(x) + /o(x) where /e(x) (^resp. fo{^)) consists of all terms of even 
(resp. odd) degree. Then 

f eH{C) if and only if fe ^ fo- 

Definition 4. If /(x) is a polynomial then /^(x) is the sum of all terms 
with highest total degree. 

Fact 3. Suppose that /(x) G HdiC) is an stable polynomial of degree n. 

(1) /^(x) is a stable homogeneous polynomial. 

(2) f^ is the limit of stable homogeneous polynomials such that all mono- 
mials of degree n have non-zero coefficient. 

(3) All the coefficients of f^ have the same argument. 
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2. Real coefficients 
Definition 5. Tid is the subset of 7id{C) whose coefficients are all positive. 

The restriction that the coefficients are all positive is natural: 

Fact 4. If fix) G H(C) has all real coefficients then all coefficients have the 
same sign. 

Proof. Suppose /(x) has degree n. We can write /(x) as a limit of /e(x) 
where /g is stable, and all coefficients of terms of degree at most n are non- 
zero. Since the signs of a stable polynomial in one variable are all the same, 
an easy induction shows that all coefficients of fe have the same sign. Taking 
limits finishes the proof. □ 

Fact 5. // a, b, c, d are positive then a + bx + cy + dxy G ^2- 

Proof. If fix, y) = a + bx + cy + dxy = where x = r + %s then 

(ac + bcT + adr + hdr"^ + bds^) 



[c + drf +d'^s'^ 



lir + ts & RHP then r > so RHP. Thus f{x, y) can not have both 
x,y in the right half plane, so / is stable. □ 

There are constructions of polynomials in 7i and 'H(C) using determi- 
nants, and they follow from the corresponding results for U. 

Fact 6. Suppose that A is skew- symmetric, S is symmetric, and all Di are 
positive definite. 

|I + XI Di + • • • + XdDd + iS\ € HdiC) 
\l + xDi + ---+XdDd + A\end 

Next we introduce polynomials that share the properties of both stable 
and upper polynomials. 

Definition 6. Let P^ = Hd H U,^. Interlacing in P^ is defined as usual 

f^giSf + ygeP+ 
H U 

and is equivalent to / < — g and / < — g 

In one variable, P^ consists of all polynomials with all real roots, no pos- 
itive roots, and all positive coefficients. Here's an example of a polynomial 
with positive coefficients that is not in P2 ■ Let /(x, y) = x{x + 3) + y (x + 
l)(x + 2). The roots of /(x, 1 + %) lie in the second and third quadrants, so 
f{x,\ + i) ^ U2 (C), and hence is not in P^. 
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3. Analytic closure 
Definition 7. 7^d(C) is the uniform closure on compact subsets of Hd{C). 
Fact 7. 

(1) y G H2d(C) and e'''^ G Wd(C). 

(2) ///(x,y) G W2d(C) t/ien e^-^^/(x,y) G W2d(C). 

(3) ///(x) G Wd(C) and ^(x) G ^^(C) i/ien f{dM^) G ^^(C) U {0}. 
This also holds for f G TidiC). 

(4) Suppose T is a non-trivial linear transformation defined on poly- 
nomials in d variables. T{e^'y) G 'H2d{C) if and only if T maps 
T-idiC) U {0} to itself 

Example 1. The Hadamard product of f{y) and g{x,y) is given by 

f *9 = i^aiy') * (^5i(x)y*) = 9i{^) 

(1) It is easy to verify that f{x) G then f{xy) G 7{2{C). 

(2) The hnear transformation exp: x* i-^ x^/il determines a map 
TCi{C) — > 7^1 (C). This follows from the above since the generating 
function f{x,y) of exp is g{xy) where g = G Pj^. 

(3) The linear transformation g ^ f * g has generating function 

= e-"exp/(yz;) 

This is in 7^(C), so we have a map x 7id{C) — > Hd{C). It is 
surprising [8j that the Hadamard also maps TCi x Tii — > Hi. 

Fact 8. Suppose that /(x,y) is a polynomial, and define T[g) = f{x,dy)g. 
The following are equivalent. 

(1) T:ndic)^nd{c)uo. 

(2) /(x,y) G7^2d(C). 

4. Positive interlacing 

If / < — g then f + (jg is stable for all a in the right half plane. If we 
restrict ourselves to just the real numbers in the right half plane we get 

positive interlacing. 

Definition 8. 

/(x) ~ (7(x) if and only if /(x) + r g{x) G P^ for all r > 

/(x) ~ 5(x) if and only if /(x) + r g{x) G Ti for all r > 

HP P 
Clearly both ~ and ~ are reflexive. In one variable / ~ 5 is equivalent 

P , 
to common interlacing. That is, / ~ 17 if and only if there is an /i G Pi so 
P P 

that h < — / and h < — g. In general it follows from Fact [5] that 
Fact 9. 



ASPECTS OF STABLE POLYNOMIALS 



5 



(1) If f,g G and there is an h such that h < — / and h < — g then 
J ~ 9- 

(2) IffeP+ thend,,{f)^d,^{f). 

H P 

Here are some elementary properties of ~ and ~. 

Fact 10. Suppose r, s are positive. 

(1) f ^ a iff rf sg. (2) f ~ g iff f + rg ~ g 

(3) If hen then fS^gifffhS. hg. U) f " g iff 9 ~ /■ 
(5) f '-^ g ^ f ^ g (6) f ~ g ^ f ~ g 
(V f^9 ^ f^9 (8) f^g ^ f^g 
(9)f^g^f"g (10) f"g ^ f^g 

Proof. The only ones that need any proof are the implication failures. 

/ ~ 9 ^ / ^ 9- If / = a^^, 5 = 1 and t > then f + tg = + t e 
H, but it is easy to check that + y ^ P^. 

/ ~ 5 =A / ■* — 9' The above example also shows this, since x +y ^ 
p p 

f ^ 9 =7^ / ■* — 9- If we let / = a;(a; + 3) and g = {x + l){x + 2) then 
f ^ g. However, we have seen that f + yg ^ P2 ■ 

□ 

Fact 11. 

(1) // E fi{^)y' e n then A(x) S /,,+2(x) /or A; = 0, 1, ... . 

(2) If f{^) + ... + yzg{^) + ...eH thenf^g. 

(3) If f < — g then f ~ xg. 

(4) If f < — g and fi < — gi then ffi ~ ggi . 

Proof. Since + i € for positive t, we know dy+t preserves TL. We first 
differentiate k times, yielding a polynomial in T-L. The constant term of 

{dl + i)[fc!/fe(x) + (fc + l)!A+i(x)2/ + (l/2)(fc + 2)!/fe+2(x)y2 + • • •] 

= [tA;!A,(x) + (fc + 2)!/,+2(x)] + y(- • • ) 

is in TL, which proves the first part. For the second part we multiply by 
yz + t where t > and consider the coefficient of yz. 

The third part follows from first part and the expansion 

(/ + y9){xy + l) = f + y{g + xf) + xgy"^ 

Next, we compute 

(/ + 2/5) (/i + y9i) = ffi + y{f9i + fi9) + y'^991 

and the conclusion follows as above. □ 
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5. Properties of one variable 

In this section we assume d = 1. When we restrict ourselves to one 
variable there are some useful techniques for showing interlacing holds. We 
let Qi denote the i'th quadrant. 



Fact 12. 










(1) / 


H 


9 




— C\(- 


(2) 


H 


- 9 




RHP. 


(3) f< 
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- 9 




-^Qi- 


(4) / 


p 


9 




RHP 


(5) / 


p 


9 


implies f 


H 

< — 9- 



Proof, f g iS f + tg e Til ioT t > iS -^(cr) (0, oo) for a e RHP iff 

i : RHP — >C\ (-00, 0) iff ^ : Qi — ^ C \ (-oo, 0) since / and g have real 
coefficients. 

Next, / ^ giSf+yg G H2iff-^(cj) RHP for x G RHP iff ^: RHP — > 
RHP iff ^ : Qi — > RHP. 

If / < — g then we take / to be monic, and we can write / = JK^ + ^i) 
and g = "^Oif /{x + rj) where n and Oj are non-negative. If cr G Qi then 
{9/f){(^) = + e Qi, and so {f/g){a) G Qi. 

li f-\-tg G P^ for all positive t then / and g have a common inter lacing [6j, 
so there is an h satisfying h < — / and h < — g. If cr G Qi then {f /g){(y) = 
By the first part {h/ f)[a) G Qi and {g/h)[a) G Q4, so {f/g){a) G 

RHP. 



The last one follows from the previous ones. 



□ 



Fact 13. 



!)///- h,- --J^fn then / ~ /i + • • • + /n. 

2) If f < — fi and g i — gi for 1 < i < n then fg -^Y. fi9i- 

f 9 ' 
f 9' 



(3) Iff^gthen 
p 



(4) /// 
stable. 



is stable. 

g then the Bezout polynomial B{x,y) 



1 

x-y 



fi^) fiv) 
9(x) 9(y) 



IS 



Proof. The first one follows from Fact [T2lf4) since -f^^''^^^"- 
The second follows from 



/ + 



^E/^^^ = Et-^uhp 

fgA^ ^ f g 
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by Fact [121 3). For the next one, write / = IK^ + ^^'^ 9 — 
where the Oj are non-negative then 

(/'9-/5')(^) = /'E 



X + T-i 



a,; 



.^2 



and this is in the lower half plane, so f'g—g'h. In particular, f'g—g'h 

is stable. 

Finally, we show that f{x)f{y) ~ B{x, y) by using the identity 

fix) fiv) 



X -riy -ri 



□ 



Example 2. Assume that fo, fi, f2, ■ ■ ■ is an orthogonal polynomial se- 
quence. The ChristofFel-Darboux formula states that 



k 



71+ 1 



fn{x) fn+l{x) 



for certain constants A;„. It follows from Fact 1131 that 

foixf + h{xf + • • • + Uxf is stable. 

We can generalize the third part of Fact [13] in several ways. The proofs 
are similar to Fact [131 - write the desired determinant in terms of the pa- 
rameters and simple geometry implies it satisfies the conditions for positive 
interlacing. 

Fact 14. 

(1) IfT.h{x)y' and0<a<2 then [f^f S 



/i afo 
f2 h 



(2) IfZhW^P^- then 



is stable. 



fk fk+1 
fk + 1 fk + 2 

(3) Suppose that Di, D2, are positive definite matrices, where Di is 
diagonal, and D2, D3 have all positive entries. If 



\ld + xDi + yD2 + zDs \ = f{x) + y g{x) + z h{x) + yz k{x) + ■ ' 

6. Questions 



then I { ^ I is stable. 



Question 1. If f ^ g or f ^ g then do f and g have a common interlacing? 
This is true for f ^ g where f,g(z . 

Question 2. If fij{x)y^z^ S P^ and r is a positive integer then is the 
following determinant stable? 

foo fio ■ ■ ■ frO 
foi fn ■ ■ ■ frl 



for fir 



fr 
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